Preprint Numbers: ANL-PHY-9217-TH-98 
lrmcTTh/9812063 

Survey of heavy-meson observables 



M. A. Ivanov*, Yu. L. Kalinovskyt and C. D. Roberts* 

* Bogoliubov Laboratory of Theoretical Physics, 
Joint Institute for Nuclear Research, 141980 Dubna, Russia 

^Laboratory of Computing Techniques and Automation, 
Joint Institute for Nuclear Research, 141980 Dubna, Russia 

^Physics Division, Bldg. 203, Argonne National Laboratory, Argonne IL 60439-4843 
(Pacs Numbers: 13.20.-v, 13.20.Fc, 13.20.He, 24.85. +p) 



Abstract 



We employ a Dyson-Schwinger equation model to effect a unified and uni- 
formly accurate description of light- and heavy-meson observables, which 
we characterise by heavy- meson leptonic decays, semileptonic heavy-to-heavy 
and heavy-to-light transitions - B — * D* , D, p, tt; D — » K* , K, it, radiative 
and strong decays - B* s ^ — > B^j; D* s ^ — > D^j, Dir, and the rare B — ► if* 7 
flavour-changing neutral-current process. We elucidate the heavy-quark limit 
of these processes and, using a model-independent mass formula valid for all 
nonsinglet pseudoscalar mesons, demonstrate that their mass rises linearly 
with the mass of their heaviest constituent. In our numerical calculations we 
eschew a heavy-quark expansion and rely instead on the observation that the 
dressed c, 6-quark mass functions are well approximated by a constant, inter- 
preted as their constituent-mass: we find M c = 1.32 GeV and Mj, = 4.65 GeV. 
The calculated heavy-meson leptonic decay constants and transition form fac- 
tors are a necessary element in the experimental determination of CKM matrix 
elements. The results also show that this framework, as employed hitherto, 
is well able to describe vector meson polarisation observables. 
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I. INTRODUCTION 



Mesons are the simplest bound states in QCD, and their non-hadronic electroweak inter- 
actions provide an important tool for exploring their structure and elucidating the nonpertur- 
bative, long-distance behaviour of the strong interaction. That elucidation is accomplished 
most effectively by applying a single framework to a broad range of observables, and in 
this the bound state phenomenology [|],|2| based on Dyson-Schwinger equations (DSEs) 
has been successful; e.g., with its simultaneous application to phenomena as diverse as tt-tt 
scattering {||| , the electromagnetic form factors of light pseudoscalar mesons , anoma- 
lous pion H and photo-pion [^jHTlJ processes, and the diffractive electroproduction [[11]] and 



electromagnetic form factors |T2[ of vector mesons. Herein we extend this application and 
implement a simplification valid for heavy quarks, so obtaining in addition a description of 
heavy-meson observables. We illustrate that by reporting the simultaneous calculation of 
a range of light-meson observables and heavy-meson leptonic decays, semileptonic heavy- 
to- heavy and heavy-to-light transitions - B — > D*, D, p, 7r; D — > K*, K, 7T, radiative and 
strong decays - BTs — ► B^j] D? s ^ — ► D( s )7, Dtt, and the rare B — > K*j flavour- changing 
neutral-current process. This is an extensive but not exhaustive range of applications. 

To introduce the heavy-quark simplification we observe that mesons, whether heavy or 
light, are bound states of a dressed-quark and -antiquark, where the dressing is described 
by the quark Dyson-Schwinger equation (DSE) ||:[] 

SAP)- 1 ■= -pA f (p 2 ) + B f (p 2 ) = A f (p 2 ) (i 7 • p + M f (p 2 )) (1) 

= ■ p + m h f ™) + Z 1 / A ^ p (p- g )y T ^( g )r^(g,p). (2) 



Here / (= u,d,s,c,b) is a flavour label, D flll (k) is the dressed-gluon propagator, Tl a (q,p) 
is the dressed- quark-gluon vertex, mj m is the A-dependent current-quark bare mass and 
J g A := J A d 4 q/ (27r) 4 represents mnemonically a translationally-invariant regularisation of the 
integral, with A the regularisation mass-scale. The renormalisation constants for the quark- 
gluon-vertex and quark wave function, Zi(( 2 ,A 2 ) and Z 2 (C 2 ,A 2 ), depend on the renormali- 
sation point, (, and the regularisation mass-scale, as does the mass renormalisation constant 
Z m (( 2 ,A 2 ) := ^(C 2 ? ^ 2 )~ 1 Za(C 2 , A 2 ). However, one can choose the renormalisation scheme 
such that they are flavour-independent. 

This equation has been much studied and the qualitative features of its solution eluci- 
dated. In QCD the chiral limit is defined by m = 0, where rh is the renormalisation-point- 
independent current-quark mass. It follows that in this case there is no scalar, mass-like 
divergence in the perturbative evaluation of the quark self energy. Hence, for p 2 > 20 GeV 2 
the solution of Eq. (0) for the chiral-limit quark mass-function is [13| 



2 ^ 2^ (-(»)") 



1 We use a Euclidean formulation with {7^,7^} = 25^, 7^ = 7^ and p ■ q = J2i=iPiQi- A vector, 



kn, is timelike if k 2 < 0. 
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where 7 m = 12/(33 — 2iV/) is the gauge-independent mass anomalous dimension and (qq)° 
is the renormalisation-point-independent vacuum quark condensate. The existence of dy- 
namical chiral symmetry breaking (DCSB) means that (qq)° 7^ 0, however, its actual value 
depends on the long-range behaviour of D^ik) and T° a (g,p), which is modelled in contem- 
porary DSE studies. (qq)° ~ — (0.24GeV) 3 is consistent with light-meson observables | fL4[ . 
In contrast, for rhf 7^ 0, 

MAP 2 ) lafg = P2 7V, p.™/, TV^T - (4) 



(I In [p 2 /A 



2 

QCD 



An obvious qualitative difference is that, relative to Eq. (|3j), the chiral-limit solution is 
l/p 2 -suppressed in the ultraviolet. 

There is some quantitative mo del- dependence in the momentum-evolution of the mass- 
function into the infrared. However, for any forms of D fJLV {k) and ^l a {q,p) that provide an 
accurate description of f^K and m nj K, one obtains Jl5| quark mass-functions with profiles 
like those illustrated in Fig. IB. The evolution to coincidence between the chiral-limit and 
u, <i-quark mass functions, apparent in this figure, makes clear the transition from the per- 
turbative to the nonperturbative domain. The chiral limit mass-function is nonzero only 
because of the nonperturbative DCSB mechanism whereas the u, <i-quark mass function is 
purely perturbative at p 2 > 20 GeV 2 , where Eq. (Q) is accurate. The DCSB mechanism thus 
has a significant effect on the propagation characteristics of u, d, s-quarks. 

However, as evident in the figure, that is not the case for the 6-quark. Its large current- 
quark mass almost entirely suppresses momentum-dependent dressing, so that M b (p 2 ) is 
nearly constant on a substantial domain. This is true to a lesser extent for the c-quark. 

We employ Cf := M^/m'j as a single, quantitative measure of the importance of the 
DCSB mechanism; i.e., nonperturbative effects, in modifying the propagation characteristics 
of a given quark flavour. In this particular illustration it takes the values: 



/ 



u,d s c b , . 

150 10 2.2 1.2 ' ^ ' 



These values are representative: for light-quarks C q = u ,d,s ~ 10-100, while for heavy-quarks 
£q =C) 6 ~ 1, and highlight the existence of a mass-scale characteristic of DCSB: M x . The 
propagation characteristics of a flavour with rrif < M x are significantly altered by the DCSB 
mechanism, while for flavours with ^> M x momentum-dependent dressing is almost 
irrelevant. It is apparent and unsurprising that M x ~ 0.2 GeV ~ Aqcd- As a consequence 
we anticipate that the propagation of c, 6-quarks can be described well by replacing their 
mass-functions with a constant; i.e., writing 

S Q (p) = - —x-,Q = c,b, (6) 

Z7 • p + M Q 

where Mq is a constituent-heavy-quark mass parameter .0 We expect that a good description 



2 Although not illustrated explicitly, when Mj(p 2 ) ~ const., Af(p 2 ) ~ 1 in Eq. (0). Eq. (Q) is an 
implicit assumption in the formulation of Bethe-Salpeter equation models of heavy-mesons, such 
as Ref. || . 
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of observable phenomena will require Mq « Mq . 

When considering a meson with a heavy-quark constituent one can proceed further, 
as in heavy-quark effective theory (HQET) |17|], allow the heaviest quark to carry all the 
heavy-meson momentum: 

P M =: m HVfl =: (Mq + E H )v„ t (7) 

and write 



H 



where k is the momentum of the lighter constituent. In the calculation of observables, the 
meson's Bethe-Salpeter amplitude will limit the range of \k\ so that Eq. @ will only be 
a good approximation if both the momentum-space width of the amplitude, lojj, and the 
binding energy, E^, are significantly less than Mq. 



In Ref. |L8| the propagation of c- and 6-quarks was described by Eq. (J8|), with a goal 
of exploring the fidelity of that idealisation. It was found to allow for a uniformly good 
description of .By-meson leptonic and semileptonic decays with heavy- and light-pseudoscalar 
final states. In that study, corrected as described below, ojb s ~ 1-3 GeV and E Bf « 0.70 GeV, 
both of which are small compared with M « 4.6 GeV in Fig. [TL so that the accuracy of the 
approximation could be anticipated. It is reasonable to expect that uj^ ujb and Ed ~ Eb, 
since they must be the same in the limit of exact heavy-quark symmetry. Hence in processes 
involving the weak decay of a c-quark (M c w 1.3 GeV) where a Dj-meson is the heaviest 



participant, Eq. (g) must be inadequate; an expectation verified in Ref. 

The failure of Eq. (|8|) for the c-quark complicates or precludes the development of a 
unified understanding of Df- and .By-meson observables using such contemporary theoret- 
ical tools as HQET and light cone sum rules (LCSRs) |T§|] . However, the constituent-like 
dressed-heavy-quark propagator of Eq. (|) can still be used to effect a unified and accurate 
simplification of the study of these observables. Herein, to demonstrate this, we extend 
Refs. [0,0 an d employ Eq. (|6|), with Mq treated as free parameters, and parametrisa- 
tions |6|-P , ITj 12[ of the dressed-light-quark propagators and meson Bethe-Salpeter ampli- 
tudes in the calculation of a wide range of observables, determining the parameters in a 
X 2 -fit to a subset of them. It is an efficacious strategy. 

Our article is divided into eight sections with a single appendix. We discuss heavy- and 
light-meson leptonic decays in Sec. [II], and their masses in Sec. |T|. In Sec. [IV] we introduce 
the impulse approximation to the semileptonic decays of heavy-mesons, and describe the 
light-quark propagators and meson Bethe-Salpeter amplitudes necessary for their evaluation. 
The impulse approximation to the other processes is presented in Sec. [V], while in Sec. |VT| we 
elucidate the heavy-quark symmetry limits of all the decays and transitions. The accuracy 
of these heavy-quark symmetry predictions is discussed in conjunction with the complete 
presentation of our results in Sec. |V1J| and Sec. |V111| contains some concluding remarks. 



II. LEPTONIC DECAYS 
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A. Pseudoscalar Mesons 



The leptonic decay of a pseudoscalar meson, P(p), is described by the matrix element [13 
f P p, := (0\Q(T p ) T j^Q\P(p)) =trZ 2 f k (T^) T 757 , XP (k; p), (9) 

J k 

where Xp(k;p) = S(k + p)T P (k;p) S(k), Q = column(u, d, s, c, b), T p is a flavour 
matrix identifying the meson; e.g., T w+ = |(A 1 + iA 2 ), (T P ) T is its transpose, S = 
diag(S u , Sd, S s , S c , Sb), and the trace is over colour, Dirac and flavour indices. T P is the 
meson's Bethe-Salpeter amplitude, which is normalised canonically according to: 

2p, = tr [ A T P (q; -p) ^±A V p{q . p) S{q ) 

Jq Op^ 
r A /-A r)K rs (a h- n) 

+ \xp{r,-p)U * [xp(k;p)U, (10) 

Jq Jk Op^ 

where: f ]b(k, — p) = C~ 1 T P (—k,—p)C, with C = 7274, the charge conjugation matrix; 
r,s,t,u are colour-, Dirac- and flavour-matrix indices; and K is the quark-antiquark scat- 
tering kernel. Equation (§) is exact in QCD: the A-dependence of Z 2 ensures that the 
right-hand-side (r.h.s.) is finite as A — > 00, and its (- and gauge- dependence is just that 
necessary to compensate that of Xp(k;p). 



The leptonic decay constants of light pseudoscalar mesons, n, K, are known [50 
f n = 0.131 GeV, fx = 0.160 GeV. The increase with increasing current-quark mass is eas- 
ily reproduced in DSE studies JTBJ and continues until at least 3m s J2TJ, at which point 



the renormalisation-group-improved and confining ladder-like truncation of K used in those 
studies becomes inadequate, and that model no longer allows the mp-dependence of f P to 
be tracked directly. However, we note from Eq. (0) that 

g P (k;p):=—r- — T P (k;p) < 00 ,m P — ► 00 ; (11) 
y/m P 

i.e., that Gp(k;p) is mass-independent in the heavy-quark symmetry limit, and hence it 
follows |22|] from the general form of the meson Bethe-Salpeter amplitude and Eqs. @ 
- ( |TT]| ) that for large pseudoscalar meson masses 

f P oc l/y/mp . (12) 

In this mo del- independent result we recover a well-known general consequence of heavy- 
quark symmetry flIT| . However, the value of the current-quark mass at which it becomes 
evident is unknown in spite of the many studies that report values of fp> and fp, some of 



which are tabulated in Ref. |23j, and the results of lattice simulations, a summary |24| of 
which reports 

f D = 200 ± 30 MeV , f B = 170 ± 35 MeV . (13) 

In Fig. we illustrate the behaviour of f P we anticipate based on these observations. It 
suggests that D- mesons lie outside the domain on which Eq. ( |T2"D is manifest. 



5 



B. Vector Mesons 



The leptonic decay of a vector meson, V\(p), is described by the matrix element (e^(p) 
is the polarisation vector: e x {p) ■ p = 0) 

fvMyefa) := (0\Q(T p f^Q\V x (p)), 

=> / y M y = |trZ 2 / A (T p ) T 7 M X^(fc;p), (14) 

where X^(^iP) — + p) r^(fc;p) «S(/c) with T^(k;p) the vector meson Bethe-Salpeter 
amplitude, which is transverse: 

Pll Tl(k;p) = 0,p 2 = -M*, (15) 

and normalised according to 

dS(q +p) 



2p,= |tr/ Tl{ q --p) a -p*±^Tl{ q -p)S{ q ) 

Jq Op^ 



dKi:(q,k;p) 



+ J q J k [x v At, -p)]„ y w [xl(k; P )U, (16) 



an analogue of Eq. fllPf) . The obvious analogue of Eq. (jTlf) is true. 

Such decays are difficult to observe directly but it is possible to estimate f p and hence 
identify the natural scale of fv- In the isospin-symmetric limit the p° — ► e + e~ decay con- 
stant, g p , is obtained from the matrix element 



9 P 

so that 



M 2 

-^e» := (0|u7 M u|pS(p)> = \ (0|« 7 ^IPa (?)> (17) 



/p = v/2^- (18) 

The experimentally measured width |2{J r p o^ e + e - = 6.77 ± 0.32 keV yields g p = 5.03 ± 0.12 
and hence 

/ p = 216±5MeV. (19) 

For the vector-mesons it follows from the general form of the Bethe-Salpeter amplitude, 
and Eqs. (|§|), flT3j) and (|TB[) that for large vector- meson masses 

fv cc l/y/My , (20) 

which again is a model-independent and general consequence of heavy-quark symmetry. In 
fact, as we illustrate below, 

f P = fv oc 1/x/wq , m Q -> oo ; (21) 
i.e., observables are spin- independent in the heavy-quark symmetry limit [17]]. 
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III. PSEUDOSCALAR MESON MASSES 



Flavour nonsinglet pseudoscalar meson masses satisfy 

f P m 2 P = M%r%, M% := tr fl „ [M c {t p , (T p ) t }] , (22) 
with = diag(m£, rrv d , m^, m£, m^) and 



i r P 



trZ 4 jf (T p ) T l5X p(k;p). (23) 

The renormalisation constant, Z 4 , ensures that r P is finite as A — > oo, and its £- and gauge- 
dependence is just that necessary to ensure that the product on the r.h.s. of Eq. fl22|) is 
gauge invariant and renormalisation point independent. 

In the limit of small current-quark masses one obtains ]13|] what is commonly called 
the Gell-Mann-Oakes-Renner relation as a corollary of Eq. (p2|); i.e., mp oc rhf, rhf — > 0. 
However, it also has an important corollary in the heavy-quark symmetry limit. Using 
Eqs. (|) and (0) 

r P oc y/mp , mp — > oo , (24) 
from which Eqs. ( |T2]) and ( |2"2"D yield 

m P oc m Q , m Q — > oo . (25) 

Thus the quadratic trajectory, valid when the current-quark mass of the constituents is small, 
evolves into a linear trajectory when this mass becomes large. In all phenomenologically ef- 
ficacious DSE models the linear trajectory is manifest at twice the s-quark current-mass |15j 
so that m^ lies on the extrapolation of the straight line joining mp, and mp in the {rhf, mp)- 
plane Hj. 



IV. SEMILEPTONIC TRANSITION FORM FACTORS 



A. Pseudoscalar meson in the final state 



The pseudoscalar — > pseudoscalar transition: P\{pi) — > ^2(^2) ^ v ', where Pi represents 
either a B or D and P 2 can be a D, K or 7r, is described by the invariant amplitude 

A(P 1 ^P 2 £u) = ^V f , f £ril->y s )uM**(p 1 ,p 2 ), (26) 

where Gf = 1.166 x lCT 5 GeV -2 , Vff is the relevant element of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix, and the hadronic current is 

M^( Pl ,p 2 ) := (P 2 (p2)\fl,f\Pi(pi)) = /+(*)(Pi+a)m + (27) 

with t := — q 2 = —(pi — p 2 ) 2 - The transition form factors, f±{t), contain all the information 
about strong-interaction effects in these processes, and their accurate calculation is essential 
for a reliable determination of the CKM matrix elements from a measurement of the decay 
width (t± := {mp 1 ±mp 2 ) 2 ): 

r(p x - p 2 tv) = ^ \v f , f \ 2 JL jT*" rft | /+ (t)| 2 [( t+ - t)(t_ - t)f/ 2 (28) 
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B. Vector meson in the final state 



The pseudoscalar — > vector transition: P(p\) — > V\(p 2 )£v, with P either a 5 or D and 
V A a D*, K* or p, is described by the invariant amplitude 



A(P ^V x £v) = ^V rf !t 1 ,(l- l5 )uet(p 2 )M^( Pl ,p 2 ) , 
where the hadronic tensor involves four scalar functions 



(29) 



,P2) 



el (m P + My) A 1 (t) + ( Pl + p 2 ) ll e x ■ q 

a Mt) A 

+ -q —TT + e„ p la p 2p 

m P + My 



Mt) 



m P + My 

2V(t) 
m P + M v 



The contribution of A 3 can be neglected unless £ = t. 
Introducing three helicity amplitudes 

H ± ^( m p + My) Mt ) T ^^lv(t), 



H = 
1 

2My^t 



m P -M v -t 



mp + My 



[mp + My] A^t) 



X(m 2 p,M v ,t) 
m P + My 



A 2 (t) 



(30) 



(31) 



(32) 



where X(mp, M v ,t) = [t+—t] [t- — t], t± = (mp±My) 2 , the transition rates can be expressed 
as 



dT 



±,o 



dt 



1927i 3 m 3 Pl 



\V f , f \ 2 t\Hm P: M v ,t)\H ±fi (t)\\ 



(33) 



and the transverse and longitudinal rates and widths are 

gTY dT + dT_ 



dt dt dt 



dt ' 



dt 



dt ' 



r-dt dTL 

Jo 



dt ' 



(34) 
(35) 



with the total width: T = Tt + Tl- The polarisation ratio and forward-backward asymmetry 
are 



a = 2 rY -1 ' Afb = 1^— 



(36) 
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C. Impulse Approximation 



We employ the impulse approximation in calculating the hadronic contribution to these 
invariant amplitudes: 



K x ( PllP2 ) = 

2N c ti D f A V x {k-,-p 2 )S q {k 2 )iOf{k 2 M)S Q {k 1 )Vp{k-, Pl )S ql {k), (37) 

J k 

where the flavour structure has been made explicit, k\p = k + Pi, 2 and: 

K x=p z( Pl ,P2) = M™( Pl ,p 2 ) , (38) 
K X=V \ Pl ,P2) = el(p 2 ) M^(p u p 2 ) ; (39) 

r x=v \(k;p) = e x (p) ■ T v (k;p)] and 0^(k 2 , k±) is the dressed-quark- W-boson vertex, which 
in weak decays of heavy-quarks is well approximated [0,0 by 

0*?(*2,*i) = 7m(1-76)- (40) 

The impulse approximation has been used widely and efficaciously in phenomenologi- 
cal DSE studies; e.g., Refs. f6Hl^, [l8l , FI . It is self-consistent only if the quark- ant iquark 
scattering kernel is independent of the total momentum. Corrections can be incorporated 
systematically and their effect in processes such as those considered herein has been 
estimated |Tj],^8|]: they vanish with increasing spacelike momentum transfer and contribute 
< 15% even at the extreme kinematic limit, t = t_. 



D. Quark Propagators 

To evaluate H PlX (pi,p 2 ) a specific form for the dressed-quark propagators is required. 
As argued in Sec. |, Eq. (|J) provides a good approximation for the heavier quarks, Q = c,b, 
and we use that herein with Mq treated as free parameters. For the light-quark propagators: 

s f (p) = -ii-v4tf) + = i7 . pA J } + ■ 

f = u,s (isospin symmetry is assumed), we use the algebraic forms introduced in Ref. ||, 
which efficiently characterise the essential and robust elements of the solution of Eq. (0) and 
have been used efficaciously in Refs. ||6|^ |lli[l2i|l8|,^2| : 



a f s (x) = 2m s T(2(x + m))) + ^(Wx^hx) (b f + b f 2 T{ex)) , (42) 
= ~ Sz I ^ 2 , (43) 



J-{y) = (1 — e y )/y, x = p 2 /\ 2 ; rrif = rrif/\; and 

4(x)=\a f s (p 2 ), 4{x) = A 2 a f v {p 2 ) , (44) 
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with A a mass scale. This algebraic form combines the effects of confinement^ and DCSB 
with free-particle behaviour at large, spacelike p 2 .[] 
The chiral limit vacuum quark condensate is [O 



(qq)° c := lim iV c tr D A 2 



S (k), 



(45) 



where at one-loop order Z 4 (( 2 ,A 2 ) = [a(A 2 )/a(( 2 )] 7m ( 1+? / 3 ), with £ the covariant-gauge 
parameter (£ = specifies Landau gauge). The ^-dependence of Z 4 (^ 2 ,A 2 ) is just that 
required to ensure that (qq)® is gauge independent. The parametrisation of Eq. (fE| ) provides 
a model that corresponds to the replacement 7 m — > 1 in Landau gauge, in which case, with 
5*0 := S u ,m=o, Eq. (||) yields 



{uu)c 



A 



4tt 2 b\b% 



In 



c 2 



A 2 



(46) 



This is a signature of DCSB in the model. 

In Ref. H the parameters %, ^ in Eqs. 



and (f43|) take the values 





rrif 


b f 


b{ 


b f 2 




it 


0.00897 


0.131 


2.90 


0.603 


0.185 


s 


0.224 


0.105 


2.90 


0.740 


0.185 



(47) 



with A = 0.566 GeV, which were determined @ in a least-squares fit to a range of light- 
hadron observables, and we note that with A QCD = 0.2 GeV they yield (ww) 1GeV 2 = 
(—0.22 GeV) 3 and (ss) 1GeV 2 = 0.8 (ww) 1GeV 2. Herein we reconsider this parametrisation 
and allow m U)S , bi' s and b\ ' s to vary. This is a reasonable step provided that in re- fitting to 
an increased sample of observables the light-quark propagators are pointwise little changed. 



E. Bethe-Salpeter amplitudes 

1. Light Pseudoscalar Mesons 

The light-meson Bethe-Salpeter amplitudes can be determined reliably by solving the 
Bethe-Salpeter equation (BSE) in a truncation consistent with that employed in the quark 
DSE [O. However, since for the light-quarks we have parametrised the solution of the quark 



3 The representation of S(p) as an entire function is motivated by the algebraic solutions of Eq. ^) 
in Refs. [25] and the concomitant absence of a Lehmann representation is a sufficient condition for 



confinement Pj3|,p6|. 

At large-p 2 : o~v(p 2 ) ~ 1/p 2 and o~s{p 2 ) ~ m/p 2 . Therefore the parametrisation does not incor- 
porate the additional lnp 2 -suppression characteristic of QCD. It is a useful simplification, which 
introduces model artefacts that are easily identified and accounted for. e = 10 -4 is introduced 
only to decouple the large- and intermediate-p 2 domains. 
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DSE, we follow Refs. |5|^ JTT|JT^J18|j2^| and do the same for the light-meson amplitude; i.e., 
for the Ti- and i\~-mesons we employ T njK (k; P) = i^S^^ik 2 ) with, 

S P (k 2 ) = ^B P (k 2 ), P = 7i,K, (48) 
Jp 

fn = /ttA/^j where Bp := B u \ bu _+ b p, obtained from Eq. (|41~D , and b^ K are allowed to vary.f] 
This Ansatz follows from the constraints imposed by the axial-vector Ward-Takahashi iden- 
tity and, together with Eqs. (f42|) and (|^), provides an algebraic representation of Xp(k;p) 
valid for small to intermediate meson energy [§,0. 

With this representation, Eqs. (^2|) and (p3|) yield the following expression for the ir- and 
.fT-meson masses: 

f 2 m 2 P = -(m u + m fP )(qq) 1 P GcV2 , (49) 
where = m^, mfK = m s , and the "in- meson condensate" is 

ml ° cV '- x3ln ^& P ^- K ' (50) 

and takes typical values |H| (qq)l Gey2 = 1.05 (uu) lGeY 2 and {qq)]^ eN = 1-64 (uu) 1GeY 2. 



2. Light Vector Mesons 



The application of DSE-based phenomenology to processes involving vector mesons is less 
extensive than that involving pseudoscalars. Therefore the modelling of vector meson Bethe- 
Salpeter amplitudes is less sophisticated. Solutions of a mutually consistent truncation 

^ indicate that a given vector meson is 



of the quark DSE and meson BSE; e.g., Ref. 
narrower in momentum space than its pseudoscalar partner but that for both vector and 
pseudoscalar mesons this width increases with the total current-mass of the constituents. 
These observations are confirmed in calculations of vector meson electroproduction cross 



sections ]TTJ and electromagnetic form factors jl2 |. The simple Ansatz 



7m +Pm 



1-V 

M 2 V/ 



(51) 



where <p{k 2 ) = 1/(1 + k 4 /ujy) with uoy a parameter and M v fixed by Eq. (|16"D, allows for the 
realisation of these qualitative features and, from Ref. [Ol, we expect lok* 



1.6 uj p . 



5 In the following we explicitly account for the flavour structure in the hadronic tensors. With 
Eq. (H) we correct an error in Eq. (33) of Ref. g§, which led to < 10% underestimates of /£*(()), 
f+ K (0) and f+ n (0). A corrected Table I, accounting also for a factor of v2 arising through a 
mismatch between the normalisation conventions for light- and heavy-mesons, is obtained with 
E = 0.698 GeV and A = 1.273 GeV, and yields Y?/N = 0.59 cf. 0.48 therein. 
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3. Heavy Mesons 



Renormalisation-group-improved ladder-like truncations of K employed, e.g., in 
Refs. are inadequate for heavy- mesons; one reason being that they do not yield the 



Dirac equation when the mass of one of the fermions becomes large. While an improved 
truncation valid in this regime is being sought, there is currently no satisfactory alternative 
and the ladder-like truncations have been used in spite of their inadequacy p9|,30|. Such 



studies cannot yield a complete and quantitatively reliable spectrum, however, the result 
that heavy mesons are described by an amplitude whose width behaves as described in 
Sec. [IV E2j must be qualitatively robust. We therefore use a simple model for the ampli- 
tudes that allows a representation of this feature: Eq. (|5lD for heavy vector mesons, with 
ip(k 2 ) — > ifnik 2 ), and its analogue for heavy pseudoscalar mesons 

r P (k;p) = j^t l5V > H (k 2 ), (52) 

where ipnik 2 ) = exp (—k 2 /ujj) and the normalisation is fixed by Eq. (|10D. We assume the 
widths are spin and flavour independent; i.e., ujb = ^b* = ^b s , etc., as would be the case in 
the limit of exact heavy-quark symmetry, which is a useful but not necessary simplification. 



V. OTHER DECAY PROCESSES 



A. Radiative Decays 

The radiative decays: H*(pi) — > H{p 2 ) j(k) , where H = -D( s ), -B( s ), are described by the 
invariant amplitude 



A(H* - H 7) = 4 B (pi)&{k) [e Q M^{p 1 ,p 2 )+e q M^(p 1 , P2 ) 
= e^ape^ipi) &(k)pi a kf, [e Q J Q (t) + e q J q {t)} , 



(53) 
(54) 



where t = —k 2 = —{p\ —P2) 2 — and ef is the fractional charge of the active quark in units 
of the positron charge. The sum indicates that the decay occurs via a spin-flip transition by 
either the heavy or light quark. The width is 



a,. 



A 3/2 (M^,m^,0) \e Q J Q (0) + e q J q (0) 



(55) 



In impulse approximation the hadronic tensors in Eq. ([53]) are: 
M%( Pl ,p 2 ) = 



2N c ti D jf r P (£;-p 2 )s Q (e 2 )iv^(e 2 ,e 1 )s Q (e l )r;(e;p l )s g (£) 

MUpx,p 2 ) = 



2N c tv D / Tp^-p^SQ^Tl^p^S^iTl^l + kjS^l + k) 



(56) 
(57) 



where Q = c or b, q = u or s and £i j2 = £ + p± t2 - 
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Fl(£i,£ 2 ) is the dressed-quark-photon vertex, which satisfies the vector Ward-Takahashi 
identity: 



Sf(£i 



2) , 



(58) 



a feature that ensures current conservation ||. This vertex has been much studied |3T| 



and, although its exact form remains unknown, its qualitatively robust features have been 



elucidated so that a phenomenologically efficacious Ansatz has emerged |32 

z£ A (4 ll) 7m + {h + 40„ ■ {h + 4) A A (£l, £ 2 2 ) + A B (£l £\) 

F{£\ 



Y, F {£l£\) 



\ [F{1\) + F(£j)] , A F (£l 



F(£j) 



(59) 
(60) 



where F = Af,Bf] i.e., the scalar functions in Eq. (p|). A feature of Eq. (|59"D is that the 
vertex is completely determined by the dressed-quark propagator. In Landau gauge the 
quantitative effect of modifications, such as that canvassed in Ref. p3 , is small and can 
be compensated for by small changes in the parameters that characterise a given model 
study P^| . The structure in Eq. (|5^) is only important for light-quarks because, using 
Eq. (H): A Q = 1, B Q = M Q , and hence 



(61) 



B. Strong Decays 

The process H*(pi) — > H(p 2 ) n(q), with p\ = —Mfj*, pi = ~ m2 H an d q 2 = — m 2 , is 
described by the invariant amplitude 

A(H* - Hn) = M»* H "( Pl ,p 2 ) := ( Pl ) p 2fl g H * H7r . (62) 

gH*Hn can be calculated even when the decay is kinematically forbidden, as for B*, and is 
sometimes re-expressed via 

9h*Htt '■= tJ- — 9h*Htt- (63) 

We calculate the coupling using the impulse approximation, Eq. (|A1|) , and this gives the 
width: 

r ^=l^% A8/a ^'<^- (64) 

We also consider the analogous decays of light vector mesons: p(P — k\ + k 2 ) — > 
7r(ki)7r(k 2 ) and K*(P) — > K(ki)%(k 2 ). For these processes the hadronic current can be 
written 

M^(k u k 2 ) = (h - k 2 ), f + {t) + P, f- (t) , (65) 
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which we calculate using the impulse approximation, Eq. ( |A2|) . The decay constant is 

9vp« = f + (t = My) , (66) 

in terms of which the widths are 



,2 



P 



C. Rare Flavour- Changing Neutral-Current Process 

The final decay we consider is the rare, flavour- changing neutral current process: B(pi) — > 
K*(p 2 )^(q), which proceeds predominantly |35j via the local magnetic penguin operator: 

Q 7j := ^M 6 F^so-^(l + 7 5 )6, (69) 

where F^ v is the photon's field strength tensor, and the effective interaction promoting this 
process, renormalised at a scale ( ~ M b , is 

H cS = -^V*V tb C 7 Q 7l , (70) 

with C 7 (M b ) = -0.299 and \V t * s V tb \ 2 = (0.95 ± 0.03) \V cb \ 2 . Consequently, the invariant 
amplitude describing the decay is 

A(B -> 7T 7 ) = V*V tb C 7 ^ #(q) (ft) M^ K "\ Pu p 2 ) , (71) 
with, at g 2 = 0, 

^(?)M / fr x * 7 (Pi,ft) = 

^(q)9BK* 1 \£^a/3PlaP2p + <W m l - Af|.) - (p x + £2)^1^) • (72) 

We calculate the coupling using the impulse approximation to the hadronic tensor, Eq. ( |A3| ) , 
in terms of which the decay width is 

/ /If 2 \ 3 f2 

rW* 7 = a em (M b )m 3 B (l - -f J ^M^7 2 l^| 2 ^ 7 - (73) 



VI. HEAVY-QUARK SYMMETRY LIMITS 

With algebraic representations of the dressed-quark propagators and Bethe-Salpeter am- 
plitudes the calculation of all observables is straightforward. In addition, one can obtain 
simple formulae that express the heavy-quark symmetry limits. We present them here, and 
in Sec. [VI]] gauge their accuracy and relevance through a comparison with the results of our 
complete calculations. 
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A. Leptonic Decays 



To begin, using Eq. (g), Eqs. ( |IUD and (|TBp with K assumed p-independent in order 
to effect a consistent impulse approximation, and Eq. fl52|) with its analogue for the heavy 
vector mesons, one finds [E22|l 



2 :=M 2 P =M 2 



VTIh Kf 
1 N r 



ma 47r 2 jo 



du<p 2 H {z) {a f s (z) + y/ua f v {z)] 



(74) 
(75) 



where z = u — 2Eh\/u, f labels the meson's lighter quark and in this section all dimensioned 
quantities are expressed in units of our mass-scale, A.Q This illustrates Eq. ([LI]). Similarly, 
using this result with Eq. (§) in Eqs. (|9p and (|i~4j), 



fp = fv = 

1 N r 



vW Kf 2 v / 2tt 2 Jo 
which is the promised illustration of Eqs. (|12D , (|20|) and (pl|). 



du (x/w - E H ) tp H {z) {a f s (z) + \y/ua f v (z)j , 



(76) 



B. Semileptonic Heavy — » Heavy Transitions 

The semileptonic heavy— > heavy transition form factors acquire a particularly simple 
form in the heavy-quark symmetry limit. From Eqs. (|27| ) and (|37D one obtains |22f : 



f±(t)=T ± £,( 



rrip 2 ± m Pl 



2 v /mp 2 m Pl 

£ f (w) = k) J dr ^ jf ^(^w) 2 er£(* w ) + \[^°v( z w) 



(77) 
(78) 



with «/ defined in Eq. ([74]), W = 1 + 2r(l - r)(u> - 1), z w = u - 2E H Ju/W and 



w 



m 2 Pl + mp 2 - £ 



(79) 



The canonical normalisation of the Bethe-Salpeter amplitude automatically ensures that 

Zf{w = 1) = 1 (80) 
and it follows [22| from Eq. ( |f8|) that 



6 Here and in the following we employ methods analogous to that described in the appendix of 
Ref. [18 1 to simplify the integrands that express observables. 
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Similar analysis applied to Eqs. (|3"D|) and ( pTj) yields 



re 



A 2 (t) = -A 3 (t) = V(t) = T + e/H. 



(82) 
(83) 



Equations flTT|), ( p2|) and (^) are exemplars of a general result that, in the heavy-quark 
symmetry limit, the semileptonic Hf — > H'^ transitions are described by a single, universal 
function: £,f{w) [36]. In this limit the functions 



are also constant 



R^w) := (1-t/U 
R 2 {w) := (1-t/t, 
1), independent of w. 



V(t) 
A 1 (t) ' 

A 2 (t) 
A 1 (t) 



(84) 
(85) 



C. Semileptonic Heavy — ► Light Transitions 



In this case the form factors cannot be expressed in terms of a single function when the 
heavy meson mass becomes large. However, some simplifications do occur. Adapting the 
analysis employed above, Eqs. (pTj) and (1371) yield 



/±(*) 

N c Kf 

4vr 2 



du (y/u — E H ) (Ph{z)£p{z) / dry/M Q 



Mr 



-#i ± H 2 



(86) 



where z x = z + 2Xr^u, X = (Mq/2) (1 - t/M%), z 2 = (1 - r)(z + ru), with similar 
expressions for Ai,A 2 , A 3 and V, Eqs. (|A4|) - (|A7|) , and the functions forming the integrands, 
H = H( z, Zi, z 2 ), also given in the appendix. 

When E H /m H — > and Mq — > m H — > oo then X = at the end point, t = t max , and 



one obtains the following simple scaling laws [36 



Pi 



f p 2 _ fj 

It also clear from Eqs. ([56]) and 



^ 

yP2 



A P2 - A Pl 

^2 ^3 



lmp 2 
m Pl 

lm Pl 
m P2 



- ( | A 7| ) that in this limit 

f—(tm&x) f+itmax) j ^S^max) ^^(^max) 

and hence one obtains a further idealisation from Eqs. ([88]): 



(87) 
(88) 

(89) 
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1? A? \jm n - {M > 

The form factors also exhibit a factorisable mass-dependence at t = but it is modulated 
by the properties of the light-quark propagator; e.g., 



/+(«) = 
In rriH N c Kf 
W 4^2 



/ du (l-E H /y/u) tp H (z)£p(z) \m f a s s (z) + za f v (z) 

}, (91) 



with the forms of A\,Ai,Az and V given in Eq. (|A13|) . The factorised dependence on 

1 /2 

hrniH/iriff is a common feature of our impulse approximation to all heavy— flight form 
factors and differs from the l/m#-behaviour obtained in LCSRs [HJ. It arises in the 
simplification of the multidimensional integrals because our model dressed-light-quark prop- 
agators evolve into free-particle propagators at large spacelike momenta. In QCD, where 
the propagators and Bethe-Salpeter amplitudes possess an anomalous dimension, we expect 
only that lnm# — > ln 7 m#, where 7 is calculable. 

D. Other Decay Processes 

In the heavy-quark symmetry limit Eq. ([55]) can be expressed in a particularly simple 
form because 

jQ(0) = J-, (92) 
M Q 

which suggests the following "classical" formulation of the width 

where we have introduced the quark magnetons: 

P/ 2M} n ' f J/(0) V ' 

In the heavy-quark symmetry limit, Mq = Mq and provides a measure of the quark's 
inertial mass. 

The expression for the H* — > Hit coupling also simplifies: 



9h*h-k = m H r- I / du (y/u- E H ) <£h(z) 2 £p{z) x 

V27T 2 Jo 



V27 

{cry a s + \a v + (§tt - E H y/u) [vs^'v ~ °V°s]} > (95) 
exposing a linear increase with the mass of the heavy meson. 
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VII. RESULTS 



The calculation of all observables is a straightforward numerical exercise. We simplify 
the integrands using the methods illustrated in the appendix of Ref. |18[] and the expressions 
we actually evaluate are similar to those presented therein with the important difference, 
however, that herein we use the constituent-like dressed-quark propagator of Eq. (|6|) for the 
c, 6-quarks, not Eq. (|||) . 

Our model has ten parameters plus the four quark masses, all of which are fixed via a 
X 2 -fit to the N ohs = 42 heavy- and light-meson observables in Tables | and |I[ This yields^ 





m f b{ b{ 




Uy 






u 
s 


0.00948 2.94 0.733 n 
0.210 3.18 0.858 K 


0.204 p 
0.319 K* 


0.515 D 
0.817 B 


1.81 c 
1.81 b 


1.32 
4.65 



(96) 



with x 2 /d.o.f = 1.75 and xV-^obs = 1-17. The dimensionless u,s current-quark masses 
correspond to m u = 5.4 MeV and m s = 119 MeV, and the parameters 6q 
were not varied, are given in Eq. 



and 63 ,s , which 



A general comparison of Eq. ([|7]) with Eq. (^) reveals 
the light-quark propagators to be little modified, with an average change in the parameters 
of only 3%, and with this parameter set 



Mi 



0.36 GeV, Mf = 0.49 GeV. 



(97) 



These values are little changed from those obtained with Eq. (|47|): Mf = 0.35, Mf = 
and are similar to the constituent-light-quark masses typically employed in quark 

0.55 GeV. Further, ujk*/u p = 1.59, which is 
as we anticipated in Sec. [IV E2|, and our 



0.48 
models 



411: M, 



u/d 



0.33 GeV and M 9 



identical to the value calculated from Ref 
simple parametrisation yields |42 
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r p = 0.69fm, /i p = 2.44^, (98) 
values of the p-meson charge radius and magnetic moment that compare well with those in 



Ref. pf: 0.61 fm and 2.69 e /(2M p ). 

For the heavy quarks we note that their fitted masses are consistent with the estimates 
of Ref. |20| and hence that the heavy- meson binding energy is large: 



E D := m D -M c = 0.67 GeV , E B := m B -M h = 0.70 GeV . 



(99) 



Eb calculated here is identical to the corrected heavy-meson binding energy obtained using 
the heavy-quark symmetry methods of Ref. [fL8|l , while Ed is 4% less. These values yield 
Ed/M c = 0.51 and Eb/Mi, = 0.15, which furnishes another indication that while a heavy- 
quark expansion will be accurate for the 6-quark it will provide a poor approximation for 



7 In the fitting we used @: V ub = 0.0033, V cd = 0.2205, V cs = 0.9745 and V cb = 0.039; and, 
in GeV, M p = 0.77, Mk* = 0.892 and, except in the kinematic factor A(mf, mj, t) where the 
splittings are crucial, averaged D- and i?-meson masses: mo = 1.99, mg = 5.35 (from rri£, = 1.87, 
m Ds = 1.97, M D * = 2.01, M D * = 2.11, and m B = 5.28, m Bs = 5.37, M B * = 5.32, M B * = 5.42). 
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the c-quark. This is emphasised by the value of uj d = ojb, which means that the Compton 
wavelength of the c-quark is greater than the length-scale characterising the bound state's 
extent. 

Having fixed the model's parameters, in Tables |T| - [V| we report the calculated values of 
a wide range of quantities not included in the fitting. Many articles report a calculation of 
some of these quantities and Refs. 



23.4 



can be consulted for tabulated comparisons. 
The calculated ^-dependence of the semileptonic transition form factors that are the 
hadronic manifestation of the b—>c,b^u,c—>-s and c — > d transitions is depicted in 
Figs. ^ and The form factors can be approximated by the monopole form 



h(t) 



h{0) 



t/h 



(100) 



with the dimensionless values of h(0) given in Tables |IV| and [V|, and hi, in GeV 2 , listed in 

Eq. qng). 



B 
B 
D 
D 



D, D* 

TC,p 

K, K* 
n,p 



f+ 



(4.63) 2 (5.73) 2 

(5.58) 2 -(21.5) 2 

(2.31) 2 (6.70) 2 

(2.25) 2 -(7.06) 2 



(4.64) 2 (4.61) 2 

(6.94) 2 (7.06) 2 

(3.09) 2 (2.78) 2 

(3.17) 2 (2.80) 2 



(101) 



As expected, in each of the channels the magnitude of the monopole mass is determined 
primarily by the heavy-quark mass, with the actual value reflecting small modifications due 
to the light-quark degrees of freedom and bound state structure. The behaviour of all the 



form factors is consistent with lattice simulations, where results are available f24"|. 



A. Fidelity of Heavy-Quark Symmetry 



The universal function characterising semileptonic transitions in the heavy-quark sym- 
metry limit, £(u0, can be obtained most reliably from B — > D,D* transitions, if it can be 
obtained at all. Using Eq. ( [77]) to extract it from f?~* D (t) we obtain 



1.08 



(102) 



which is a measurable deviation from Eq. (p0|). We plot ^ + {w)/^+(0) in Fig. and compare 
it with two experimental fits f3~8| : 



aw) = i- P 2 (w-i), P 

^(w) = — ^— exp [(1 - 2p 2 ) 
w + 1 L 



0.91 ±0.15±0.16 



w — 1 
w + 1 



1.53 ±0.36 ±0.14. 



(103) 
(104) 



The agreement obtained here is possible because, unlike Ref. jl8| , we did not employ a 
heavy-quark expansion for the c-quark. Our calculated result is well approximated by 



(w 



1.98. 



(105) 
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Pi, 



We have also used Eqs. ([82]) and fl33|) to extract from 5 — > D* . This yields 

£ Al (l) = 0.987, = 1.03, £ y (l) = 1.30, (106) 

an w-dependence well-described by Eq. ( |105| ) but with 

= 1.79, /^ 2 = 1.99, p^ = 2.02, (107) 

and the ratios, Eqs. fl8"4] ) and (^), 

Ri{l)/Ri{w max ) = 1.08 , R 2 {l)/R 2 (w max ) = 1.06 . (108) 

This collection of results furnishes a measure of the degree to which heavy-quark symmetry 
is respected in b —>■ c processes. Combining them it is clear that even in this case, which 
is the nearest contemporary realisation of the heavy-quark symmetry limit, corrections of 
< 30% must be expected. 

The scaling laws in Eqs. (|87|) and (0), which relate the heavy— > light form factors at 
t = t_, can be tested in B, D — > ir, p decays and we find 





Imp 
V rn B 


= 0.61 






A? 
A? 


A D A D 


v D 
v B 


A? 
Af 


A° 
Ai 


f? 

n 


Ed 
M c 
E D 
M c 


= 0.51 
= 0.06 


E B _ 

M b 

E B — 
M b 


0.15 
0.02 


0.32 
0.45 


0.77 
0.82 


1.33 
0.86 


1.29 
0.87 


0.79 
0.72 


1.68 
0.99 


0.58 
0.62 



(109) 



The first line is obtained using the actual values of our model parameters, Eq. (0), while the 
second uses artificially inflated constituent-quark masses: M c = 1.76 GeV, Mj, = 5.18 GeV, 
and M p = 0. This comparison indicates that the scaling laws fail in our model because c- 
and 6-quark masses that are consistent with contemporary estimates [[21J do not validate the 
approximations Eh / rnn = = M p j mn used in the derivation of these laws. 

As noted in Sec. |V1D| , in the heavy-quark limit the radiative decays can be used to define 
an inertial-mass for the quarks. We find 





M Q J Q (0) M^ n (GeV) M,; n (GeV) 


D* - 


-> D 


0.56 2.4 1.5 


B* - 


-> B 


0.70 6.7 1.7 


D* s - 


+ D S 


0.85 1.5 2.3 


B*- 


^B s 


1.1 4.2 2.7 



(110) 



Comparing these results with Eqs. ( p6| ) and ( p7|) it is clear and unsurprising that M m is not 
a good measure of the constituent mass when the binding energy is large. 



VIII. EPILOGUE 



We have described a direct extension of Dyson-Schwinger equation (DSE) based phe- 
nomenology to experimentally accessible heavy-meson observables. In doing this we explored 
the fidelity of a simple approximation, Eq. ([]), to the dressed-heavy-quark propagator. In 
our framework this approximation is a necessary precursor to effecting a heavy-quark ex- 
pansion. However, in contrast to Refs. [I8|,22|, herein we elected not to consummate that 
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expansion and thereby achieved a unified and uniformly accurate description of an extensive 
body of light- and heavy-meson observables. In doing so our results indicate that corrections 
to the heavy-quark symmetry limit of ^ 30% are encountered in b — > c transitions and that 
these corrections can be as large as a factor of two in c — > d transitions. 

Our calculation of the semileptonic transition form factors for B- and D-mesons on their 
entire kinematic domain and with the light-quark sector well constrained is potentially useful 
in the experimental extraction of the CKM matrix elements V c b, V u b. That is true too of 
our calculation of the leptonic decay constants; e.g., accurate knowledge of fp can assist 
in the determination of V t d- They also indicate that .D^-mesons do not lie on the heavy- 
quark 1/y^mQ-trajectory. In elucidating a mass- formula valid for all nonsinglet pseudoscalar 
mesons, we demonstrated that in the heavy-quark limit pseudoscalar meson masses grow 
linearly with the mass of their heaviest constituent; i.e., nip oc itiq. Our calculations also 
provide an estimate of the total width of the DT^- and D*°-mesons, for which currently 
there are only experimental upper-bounds. 

Although this study is a significant improvement and extension of Refs. [|l"8| .p2| 



more 



is possible. One simple step is a wider study of light vector meson observables so as to 
more tightly constrain their properties. Existing models for the Bethe-Salpeter kernel are 
applicable to these systems and studies akin to Ref. JT3] are underway. A pleasing aspect of 
our study, however, is the demonstration that DSE phenomenology as applied extensively 
hitherto is well able to describe vector meson polarisation observables. A more significant 
extension is the development of a kernel applicable to the study of heavy-meson masses. That 
would provide further insight into the structure of heavy-meson bound state amplitudes. 
They are an integral part of our calculations but only rudimentary models are currently 
available. It would also assist in constraining DSE phenomenology via a comparison with 
calculations and models of the heavy-quark potential. 
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APPENDIX A: COLLECTED FORMULAE 



The impulse approximation to the hadronic tensor describing the strong decay of a heavy 
vector meson is (£\ = £ — p±, q = Pi — P2, q 2 = — Tn%.) 



K *(pi>j*)= ( A1 ) 

2V2N c tr D ^ T P (£; -p 2 ) S Q (h) irf(l; Pl ) S u (£) f n (£; -q) S u (£ + q) , 
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and that for a light vector meson is similar: 
M^(k u k 2 ) = 

2N c tr D f A iT v Jt- P) S q (i ++ ) T P (£ 0+ ; -h) S u (£ + .) r w (^_ ; -k 2 ) S u (£^) , 



(A2) 



where P = hi + k 2 , £ a p = t + f &i + §^2, and g = w, s for 1/ = p, K*. 

The impulse approximation to the hadronic tensor describing the rare neutral current 
process is (l x 2 = £ + pi 2 ) 



(A3) 



2X c tr D / r* {t,-p 2 )S s {h)q p( T w {l + l5 )S h {h)V S {t, Pl )S q {l) 



The heavy-quark symmetry limits of the leptonic heavy — > light-meson transition form 
factors are Eq. (R6|) and 



N c Kf 

Ait 2 Jo 



du (y/u - E H ) <f H (z) <pf( z ) 







A 2 (t) 
N c k s 
2tt 2 

N c K f 



du (y/u - E H ) Vh{z) 1 (z) / dr y/M Q 



1 1 

dr - — Ha-, 

VM Q 

1 



2tt 2 Jo 



du (y/u - E H ) (f H (z) <~Pi (z) / dry/M Q 



1 

Mr 



H A . 2 - if 



A3 



V(t) 



N c K f r°° 



2tt 2 Jo 



du (y/u - E H ) <Ph(z) <fY( z ) / dr \/ M Q #2 



(A4) 

(A5) 

(A6) 
(A7) 



where: z = u- 2E H \/u, z x = z + 2Xrv/M, X = (Mq/2) (1 - t/M, 2 ,), z 2 = (1 - T )(z + ru), 
and 



#1 = CTS^S -Ty/u {a S CT V - <7 V d S ) + ZOyOy , 

H 2 = a s a v + z 2 (a s a' v - a v a' s ) + Ty/u a v a v , 
H Al = cr s vs + X (a s a v + z 2 [a s a v - a v a s ]) 

-Ty/u (a s a v - o v a s ) + (Xt^u + z + z 2 ) a v a v 
Ha 2 = -2ry/u (a s a v + Ty/uOyOy) , 
Ha 3 = H 2 + 2ry/u z 2 a v a' v , 



(A8) 
(A9) 

(A10) 
(All) 
(A12) 



with: a = cr(z), a = (z x ), a' = da(z)/dz and a' = dcr(zi) / dz\. 

At t = the behaviour of these form factors simplifies further, as described by Eq. fl9T|) 



and 



4A 1 {0) = A 2 (0) = -A 3 (0) = V(0) 



\nrriH N c Kf 



m 



H 



du (1 - E H I y/u) if H {z) y{{z) a f s (z) . 



(A13) 



In Eqs. (|A4|) - ( |A13|) all dimensioned quantities are expressed in units of our mass-scale, A. 
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FIG. 1. Quark mass function obtained as a solution of Eq. (p|) using D fiu (k) and Fl a (q,p) 
from Ref. (l3j and current-quark masses, fixed at £ = 19 GeV: iTi ud = 3.7 MeV, = 82MeV, 

p 2 define the Eu- 



0.58 GeV and m\ = 3.8 GeV. The indicated solutions of M 2 (p 2 ^ 



clidean constituent-quark mass, Mf, which takes the values: M® = 0.56 GeV, Mf = 0.70 GeV, 
Mf = 1.3 GeV, Mf = 4.6 GeV. 
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2. Experimental values of fn,Ki filled circles; lattice estimates of Jd.b in Eq. (|l3|), 
open circles; our calculated values of f-K,K,D,B-> diamonds, see Sec. VI] . (We estimate a theo- 
retical error of 10%.) The dashed line is a fit to the experimental values and lattice estimates: 



fp = (0.013 + 0.028 mp)/(l + 0.055 m P + 0.15 m 2 p ), which exhibits the large-m P limit of Eq. @ 
and the dotted line is the large-mp limit of this fit. 
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FIG. 3. Upper panel: calculated form of the semileptonic B — > D and B — > D* form factors. 
Lower panel: the semileptonic B — > tt and B — > p form factors with, for comparison, data from 
a lattice simulation |39| and a vector dominance, monopole model: ff^ n (t) = 0.46/(1 — t/m 2 B *), 
mg« = 5.325 GeV, the light, short-dashed line. Monopole fits to our calculated results are given in 



Eqs. fllOOD and ( gOg ). 
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FIG. 4. Calculated form of the semileptonic D — > and D — > X* (upper panel), and 
D — > 7r and D — > p (lower panel) form factors. Monopole fits to our calculated results are given in 
Eqs. QlOOD and (ggg). 
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FIG. 5. Calculated form of £(u>) compared with recent experimental analyses. Our result: 
solid line. Experiment: short-dashed line, linear fit from Ref. [S3], Eq. (103); long-dashed line, 



nonlinear fit from Ref. [38|, Eq. (|104j ). The two light, dotted lines are this nonlinear fit evaluated 
with the extreme values of p 2 : upper line, p 2 = 1.17 and lower line, p 2 = 1.89; data points, Ref. [46]. 
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TABLES 



TABLE I. The 16 dimension- GeV quantities used in x 2 -htting our parameters. The weighting 
error is the experimental error or 10% of the experimental value, if that is greater, since we expect 
that to be the limit of our model's accuracy. The values in the "Obs." column are taken from 



Refs. [H,20|,|4|. 



Obs. 



Calc. 



Obs. 



Calc. 



U 
Ik 

(uu) 1 / 3 

f P 

r 

pTVTT 

fn 

fB 



0.131 
0.160 
0.241 
0.245 
0.216 
0.151 

0.200 ± 0.030 
0.170 ± 0.035 



0.146 
0.178 
0.220 
0.255 
0.163 
0.118 
0.213 
0.182 



m K 
(q<i}k 3 

fK* 

Tk*(Kw) 

I'D s 

g B K*fM b 



0.138 
0.496 
0.227 
0.287 
0.244 
0.051 

0.251 ± 0.030 
2.03 ± 0.62 



0.130 
0.449 
0.199 
0.296 
0.253 
0.052 
0.234 
2.86 



TABLE II. The 26 dimensionless quantities used in x 2 -fitting our parameters. The weighting 
error is the experimental error or 10% of the experimental value, if that is greater, since we expect 
that to be the limit of accuracy of our model. The values in the "Obs." column are taken 
37-f40[, The light-meson electromagnetic form factors are calculated in impulse 
and £(w) is obtained from f+^ D (t) via Eq. (fT7|). 



from Refs. [£0],| 
approximation 





Obs. 


Calc. 




Obs. 


Calc. 


tf-"(0) 


0.73 


0.58 




0.44 ± 0.004 


0.44 


F tt (3.3 GeV 2 ) 


0.097 ± 0.019 


0.077 


B(B -» D*) 


0.0453 ± 0.0032 


0.052 


P 2 


1.53 ± 0.36 


1.84 




1.25 ± 0.26 


0.94 


ea.i) 


0.86 ± 0.03 


0.84 


A B—rD* 

FB 


0.19 ± 0.031 


0.24 


£(1.2) 


0.75 ± 0.05 


0.72 


B(B -> tt) 


(1.8 ± 0.6) xl0 - 4 


2.2 


£(1.3) 


0.66 ± 0.06 


0.63 


J + (14.9 GeV 2 ) 


0.82 ± 0.17 


0.82 


£(1-4) 


0.59 ± 0.07 


0.56 


J +(n.9GeV 2 ) 


1.19 ± 0.28 


1.00 


£(1-5) 


0.53 ± 0.08 


0.50 


• y +(20.9GeV 2 ) 


1.89 ± 0.53 


1.28 


B(B -> D) 


0.020 ± 0.007 


0.013 


B(B p) 


(2.5 ± 0.9) xlo -4 


4.8 


B(D K*) 


0.047 ± 0.004 


0.049 


f^ K (o) 


0.73 


0.61 


V(0) D ^ K * 
Ai(0) 


1.89 ± 0.25 


1.74 


f D-*« (0) 


0.73 


0.67 


V L D^K* 

a-m d ^ k 

A 1 (0) 


1.23 ± 0.13 


1.17 


9B*Bir 


23.0 ± 5.0 


23.2 


0.73 ± 0.15 


0.87 


9D*Dtt 


10.0 ± 1.3 


11.0 
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TABLE III. Calculated values of a ran 
parameters, with the light-meson electroma; 
tion H||. The "Obs." values are extracted from Refs. |^||J2|^3]. 



e of observables not included in fitting 
gnetic form factors calculated in impulse 



the model's 
approxima- 





Obs. 


Calc. 




Obs. 


Calc. 


fKTK 


0.472 ± 0.038 


0.46 


~fK r K0 


(0.19 ± 0.05) 2 


(0.10) 2 


Qp'K'K 


6.05 ± 0.02 


5.27 


T D *o (MeV) 


< 2.1 


0.020 


9K*Kn 


6.41 ± 0.06 


5.96 


T D , + (keV) 


< 131 


37.9 


9p 


5.03 ± 0.012 


5.27 


T D * Dsl (MeV) 


< 1.9 


0.001 


fn* (GeV) 




0.290 


T B *+ B+1 (keV) 




0.030 


fo t (GeV) 




0.298 


r jB ,o jB o 7 (keV) 




0.015 


f Bs (GeV) 


0.195 ± 0.035 


0.194 


rs|B s7 (keV) 




0.011 


f B * (GeV) 




0.200 


B(D* + -> D+tt°) 


0.306 ± 0.025 


0.316 


f B; (GeV) 




0.209 


B(D* + ^ D°tt+) 


0.683 ± 0.014 


0.683 


fDs 1 Id 


1.10 ± 0.06 


1.10 


B(D* + ^ 


01 1 +0 - 021 

U.Ull _ .007 


0.001 


fB s 1 !b 


1.14 ± 0.08 


1.07 


B(D*°^D°tt°) 


0.619 ± 0.029 


0.826 


Id* 1 Id 




1.36 


B(D* ^ D°~f) 


0.381 ± 0.029 


0.174 


fB* 1 Ib 




1.10 


B(B -» K*7) 


(5.7 ± 3.3) 10 -s 


11.4 



TABLE IV. Calculated values of some b — > c and b - 
included in fitting the model's parameters. The "Obs." 



-> u transition form factor observables not 
values are extracted from Refs. 2C,23|. 



Obs. Calc. 


Obs. Calc. 


Af~ 


^(0) 


0.57 


At 


Vmax / 


0.88 


Af- 


-^(0) 


0.56 


A? 


^D*( t B-+D*\ 


1.16 


V B ~ 


^(0) 


0.70 


yB- 


->D* (xB—tD* \ 
rmax / 


1.47 


Rf~ 


^•(1) 


1.30 ± 0.39 1.32 


R? 




0.64 ± 0.29 1.04 


Rf~ 


" D * (Wmax) 


1.23 


R§~ 


('"max) 


0.98 


a B - 




0.60 


ft 


*D (t B->D\ 
Vmax ) 


1.21 


ft 


-(o) 


0.45 


ft 


>TT(fB-*TT\ 

Vmax / 


3.73 


aT 


*'(0) 


0.47 


At 


-*Pr t B-*p\ 
Vmax / 


0.45 


A? 


*>(0) 


0.50 




^PuB^p\ 
Vmax ) 


0.81 


v B - 


^(0) 


0.68 


yB- 


-+p( t B-*p\ 
Vmax / 


1.17 






1.15 


R B 




0.80 






1.44 


Rf 




1.06 
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TABLE V. Calculated values of some c — > s and c - 
included in fitting the model's parameters. The "Obs." 



d transition form factor observables not 



values are extracted from Refs. [23,20]. 





Obs. 


Calc. 






Obs. 


ualc. 


D{L)^ — ► p u ) 




0.032 


D— 

a 


* () 
r 

->P°) 




1 no 

1.03 


d[JJ — ► K ) 


0.03 f ± 0.002 


0.036 


B(D 
B(D- 


0.044 ± 0.034 


0.065 




0.56 ± 0.04 


0.46 




*K* ( f D^K*\ 

\ TY1 SI Y / 
\ lllcwv / 


0.66 ± 0.05 


0.47 


^->**(0) 


0.39 ± 0.08 


0.40 




+K* ( f D^K*\ 
V^max / 


0.46 ± 0.09 


0.44 


y D ^*(o) 


1.1 ± 0.2 


0.80 




*K* (+D^>K*\ 
V-max ) 


1.4 ± 0.3 


0.92 






1.72 


R§~ 






0.83 






1.74 


R§~ 


(^max) 




0.87 


B(D°^tt) 


0.103 ± 0.039 


0.098 




K (j.D—*K\ 
\ max / 


1.31 ± 0.04 


1.11 


/+ "XO) 
/J-*(0) 


1.2 ± 0.3 


1.10 


f?~ 


7T/J.Z) — >7T\ 

v'Tnax ) 




2.18 


Af^'(0) 




0.60 


A?~ 


V max l 




0.58 






0.54 


At 


rmax / 




0.64 


v D ^p(o) 




1.22 


yD- 


+p( t D->p\ 
rmax / 




1.51 






2.08 


R 2 


; p (i) 




0.88 






2.03 


R 2 


(^max) 




0.91 
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